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We explain the observed specific heat anomaly (and hence entropy change) in the colossal mag- 
netoresistive manganite Smo.55Sro.45Mn03, by introducing phase separation into the current carrier 
density collapse theory via the notion of the ferromagnetic volume fraction. Within the same frame- 
work, we have also been able to explain the observed electrical resistivity of Smo.55Sro.45Mn03 by 
using appropriate expressions governing the scattering mechanisms far away from the transition. 
Fitting specific heat and resistivity results has allowed us to separate the hopping activation energy 
of polaronic carriers and the bipolaron binding energy contributions to the exponential behaviour 
of resistivity in the paramagnetic phase. 
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I. INTRODUCTION 

The existence of a ferromagnetic-paramagnetic and 
metal-insulator transition in lanthanum manganites was 
established in the early 1950's [U [2] and was extensively 
studied thereafter. The magnetic transition which is 
associated with unusual transport properties, including 
large negative magnetoresistance, is observed in a fam- 
ily of doped manganites of perovskite structure, with the 
chemical formula Rei_a;Aa;Mn03, where Re is rare earth 
(La, Pr, Nd, Sm...), and A is a divalent metal (Ca, Sr, 
Ba...). The recent resurgence of interest in these sys- 
tems is related to the demonstration of a very large neg- 
ative magnetoresistance in thin films [21 3] , termed colos- 
sal magnetoresistance (CMR), which shows prospects of 
possible technological applications. 

The metal-insulator transition in manganites has tra- 
ditionally been attributed to the double exchange mech- 
anism, which results in a varying bandwidth of electrons 
in the Mn^^ d shell as a function of temperature and 
doping level[5]-[7]- More recently it has been realized 
[g Ej that the effective carrier-spin exchange interaction 
of the double-exchange model is too weak to lead to a 
significant reduction of the electron bandwidth. Millis et 
al[51 proposed that there is in addition to the double ex- 
change mechanism, a strong electron-phonon interaction 
arising from the Jahn- Teller splitting of the outer Mn d 
level, that plays a crucial role in the physics of mangan- 
ites. X-ray-absorption fine structure of La2/3Cai/3Mn03 
reveals differences between the nearest-neighbour Mn-0 
bond distributions at 80 K and above T;,. The single- 
site distribution at low temperature becomes a complex, 
multi-site distribution above Tc. The observed change 
is consistent with proposed small polaron related local 
Jahn- Teller distortions for T>Tc [T3]. A giant isotope 
effect [ini [H] suggest a very large electron phonon inter- 
action, the sign anomaly of the Hall effect and the Arrhe- 
nious behaviour of the drift and hall mobilities, and the 
fact that polaron hopping satisfactorily accounts for resis- 



tivity in the paramagnetic phase |12j , verify the polaronic 
nature of the charge carriers in perovskite manganites. 

Low-temperature optical [HHIS], electron energy-loss 
(EELS) [nillHj, photoemission [Ullin], and thermoelec- 
tric [2T| measurements showed that the ferromagnetic 
phase of manganites is not a conventional metal and 
confirmed that manganites were, in fact, charge-transfer 
doped insulators. EELS[T71[TH] and O Is x-ray absorp- 
tion spectroscopy [52] consistently show oxygen p holes 
as current carriers as opposed to d Mn'^"'' electrons pro- 
posed in the double exchange theory. Moreover Hall ef- 
fect measurements by Matl et al [55] on Lai_a;Caa;Mn03 
have shown the carrier mobility near Tc to be field inde- 
pendent. Also there are known classes of materials that 
exhibit colossal magnetoresistance in which it is guaran- 
teed that the double exchange mechanism is non-existent 
such as pyrochlore manganites |24j and chromium spinels 
[25] . Like the perovskite manganites these compounds 
exhibit very large reductions in their electrical resistiv- 
ity upon application of a magnetic field. However, unlike 
the perovskites there is no possibility for mixed valency, 
a prerequisite for the double exchange mechanism [26) . 
These observations strongly suggest that the double ex- 
change mechanism is not the origin of colossal magne- 
toresistance. 

These observations led to a novel theory of the 
ferromagnetic-paramagnetic phase transition driven by 
non-degenerate polarons in doped charge-transfer mag- 
netic insulators, the so-called current-carrier density col- 
lapse (CCDC) theory [27]. The CCDC theory takes into 
account the tendency of polarons to form local bound 
pairs (bipolarons) as well as the exchange interaction of 
p polaronic holes with d electrons. These bipolarons are 
practically immobile in manganites because of the strong 
electron-phonon interaction, in contrast with the case for 
cuprates, where bipolarons are mobile and responsible for 
in-plane transport ^28j. The argument behind the CCDC 
theory is that in the paramagnetic phase a large frac- 
tion of polarons are bound into immobile bipolarons. As 
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temperature decreases within the paramagnetic phase, 
the density of bipolarons increases, resulting in fewer 
mobile polarons, and hence, the resistivity quickly in- 
creases, reaching a maximum at Tc. At some tempera- 
ture Tc non-degenerate polarons polarize the manganese 
spins via the p-d exchange interaction, Jpd, one of the 
polaron sub-bands falls bellow the energy of the bipo- 
larons, thus causing pairs to break up, and the density 
of charge carriers (mobile polarons) increases suddenly. 
The occurrence of a deep minimum in mobile polaron 
density, near Tc is the cause of the peak in resistivity, 
application of an external magnetic field will also cause 
a break up of bipolarons resulting in large negative mag- 
netoresistance. The CCDC theory can explain physics 
of all colossal magnetoresistive materials including the 
pyrochlore manganites and chromium spinels. 

Hall data by Westerburg et al [22] has indeed shown 
a current carrier density collapse. Westerburg et al 
[29] analysed the Hall data of Lao.eeCao.ssMnOa and 
Lao.eeSro.asMnOa it was shown in both samples that with 
increasing temperature there is a sharp drop in the car- 
rier density near the metal-insulator transition, further 
increase of temperature causes an increase in the den- 
sity of carriers, such an observation is well explained by 
the CCCDC theory. A similar drop in the carrier den- 
sity at the metal insulator transition has been observed 
in Ndo.sSro.sMnOs [30]. There is other independent ex- 
perimental evidence for the existence of bipolarons in 
perovskite manganites. Thermoelectric power measure- 
ments by Zhao et al [HI], show an isotope effect. If small 
polarons are bound to impurity centres, there would be 
no isotope effect in thermoelectric power measurements, 
whereas, if small polarons are bound into localized bipo- 
laron states one would observe such an isotope effect, as 
the bipolaron binding energy is dependent upon the mass 
of isotopes. Perring et al ^32^ have found short range 
antiferromagnetic order above T^ in a so-called double 
perovskite, which may be evidence for the presence of 
singlet bipolarons. 

A combination of low-energy electron diffraction and 
angle-resolved photoemission spectroscopy on some semi- 
conducting interfaces has provided unambiguous evi- 
dence of a bipolaronic insulating state [33] jH] . There is 
also evidence to support the existence of immobile spin- 
singlet bipolarons in barium titanate perovskites |351 136] 
which are of a very similar structure to the CMR per- 
ovskite manganites. Furthermore x-ray and neutron scat- 
tering measurements directly demonstrate the existence 
of short-range correlations of polarons in the paramag- 
netic phase of colossal magnetoresistive perovskite man- 
ganites [STHiO] . The polaron correlations are shown to 
have correlation lengths of approximately 10A[40]. These 
short-range polaron correlations are shown to grow with 
decreasing temperature, but disappear abruptly at the 
ferromagnetic transition [37l |38l [40]. It is established 
that the temperature dependence of the polaron corre- 
lations is intimately related to the transport properties 
of the manganites [3Sl 130] • Also the polaron correlations 



collapse under an applied magnetic field [371 SO] ■ The po- 
laron correlations are interpreted as bipolarons by Nelson 
et al. [SO]- The experimental observations of polaron cor- 
relations in the paramagnetic phase of perovskite man- 
ganites are fully consistent with the formation of bipo- 
larons within the CCDC theory. 

Chipara et al fH] claim that bipolarons are absent 
in, Lao.esCao.asMnOa [41] as analysis of ESR measure- 
ments has not shown the presence of bipolarons. However 
Chipara et al [H] also accept that low-spin bipolarons 
(singlets) have a zero spin and, accordingly, cannot be di- 
rectly inferred from ESR data. We assume, as is usually 
the case, that the triplet states of bipolarons lie higher 
in energy than the singlet state, we therefore neglect the 
case of triplet bipolarons. The argument of Chipara et 
al [4lj regarding the absence of bipolarons, therefore be- 
comes irrelevant. 

While there is somewhat of a lack of spectroscopic ev- 
idence for the presence of bipolarons in perovskite man- 
ganites, given all of the experimental evidence regarding 
the observed carrier density collapse and correlated po- 
larons that are observed only above Tc the idea of such a 
bipolaronic insulating state in perovskites does not seem 
unreasonable. We believe that the CCDC theory is a 
much more complete theory of the colossal magnetore- 
sistance phenomena, and other properties of these mate- 
rials, when compared with the double exchange model. 
In this present paper we aim to show how the micro- 
scopic CCDC theory may describe the experimental spe- 
cific heat and resistivity measurements in manganites, 
when used in conjunction with a generic phase separa- 
tion model. 



II. CCDC MODEL AND THE INTRODUCTION 
OF PHASE SEPARATION 

Alexandrov and Bratkovsky [27| have proposed that 
the physics of perovskite manganites including the colos- 
sal magnetoresistance effect, can be explained by the cur- 
rent carrier density collapse theory (CCDC). 

The grand thermodynamic potential of such a sys- 
tem has been calculated and has contributions from po- 
larons, bipolarons, localized Mn^'^ spins, and the double- 
counting term, respectively [571 131]; 
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Normalizing the grand thermodynamic potential of 
the system by using the dimensionless temperature t = 
2kBT/ JpdS, magnetic field h = 2^b^I Jpd^ and bipo- 
laron binding energy 5 ~ A/Jp^S, components of the 
grand thermodynamic potential are found as |27[ H^ : 



rip — — 2i/tycosh 



(2) 



3 



BP 



fis = -t In 



-t ln[l 



sinh[(5/2)(m + 4h)/2t] 
sinh[(l/2)(m + 4h)/2t]_ 



(3) 



(4) 



V 



-erfc 



t-tc 

r 



(9) 



Here we have taken the Mn spin, S = 2 due to strong 
Hund's couphng and the degeneracy factor, is 3 due 
to the p-character of the polarons. If the bound pairs 
are extremely local objects, i.e. two holes on the same 
oxygen, then they will form a singlet. However, if these 
holes are localized on different oxygens, then they may 
well have parallel spins and form a triplet state. The 
triplet state is separated from the singlet state by an 
exchange energy, Jgt-We shall consider the simple case in 
which the separation of the triplets from the singlets, Jst 
is much larger than the critical temperature. In this case 
D=l. 

From the grand thermodynamic potential a series of 
mean field equations describing polaron density, chemical 
potential, manganeese magnetization, and polaron mag- 
netization, n, (7 and m, respectively has been found 
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where y = e^^^ The polaron density can be found as a 
function of a and temperature t by forming a quadratic 
equation by substituting Eq( 6]) into Eq([5| . However one 
is left with a transcendental function to solve for a, to 
overcome this one may utilize the fact that the transition 
is first order in a homogeneous system in a wide range 
of 5. One can assume that the magnetization behaves as 
a step function that can become "smeared" due to dis- 
order. Alexandrov et al [iSJ suggested that one should 
consider that different parts of a particular sample will 
exhibit differing values of 5 and therefore have different 
transition temperatures, and that a Gaussian distribu- 
tion would provide a reasonable model for the spread 
of these transition temperatures. One can therefore say 
that the normalized magnetization (which we call the vol- 
ume fraction, V) is equal to 1 minus the cumulative distri- 
bution function of transition temperatures. The volume 
fraction is then given by; 



r = where s is the standard distribution of transi- 
tion temperatures across the sample. 

We term the normalized magnetization as the volume 
fraction as it describes the fraction of the sample that 
remains in the ferromagnetic phase as the temperature 
is increased. 

Broadening of the magnetic transition is in agreement 
with numerical studies by Mercaldo et al [U] who have 
shown that introducing quenched disorder, such as grain 
and crystalline boundaries into a three-dimensional pure 
system that exhibits a first order transition causes broad- 
ening of the temperature dependence of magnetization. 
Bulk samples of perovskite manganites usually exhibit 
crystallite size in the order of nanometers [iSl 148| . thus 
leading to nanoscale phase separation in manganites |47j . 



III. TRANSITION TEMPERATURE AND 
VOLUME FRACTION OF Smo ssSro 45Mn03 

Egilmez et al [48l |49] have produced samples of 
Smo.55Sro.45Mn03 annealed at difii'erent temperatures to 
produce different grain and crystallite size. Annealing at 
higher temperatures produced larger grains and crystal- 
lites, reducing the density of grain boundaries, and hence 
reduced disorder. The extra disorder in the low tempera- 
ture annealed samples results in broadening of the tran- 
sition and hence volume fraction, see Fig([2]). The erfc 
function, EqQ corresponds to a Gaussian distribution of 
transition temperatures, however the magnetization data 
of Smo.55Sro.45Mn03 can be better fitted using combina- 
tions of polynomials and the erfc function. In general one 
can use any expression to describe the volume fraction, 
as in Fig([2]). 

One should remember that within the CCDC model 
that although the density of polarons is temperature and 
field dependent the total number of holes doped into the 
system is constant and equal to the doping concentra- 
tion, N=x. Utilizing the fact that Helmholtz free energy 
is continuous over a phase transition we may calculate, 
F = f2 -|- /xN in both the ferromagnetic and paramag- 
netic phases, by substituting tr = 1 or respectively 
into the expressions for f2, and /i. The temperature 
at which these two functions cross over is the transi- 
tion temperature. FigQ shows a plot of free energy in 
the CCDC model for both phases, with 5 = 0.65 and 
x= 0.45. This value of 5 was chosen such that a reason- 
able agreement between experimental specific heat data 
of Smo.55Sro.45Mn03 [48, .49, and CCDC theory can be 
achieved (section(IV)). We have assumed 5 = 0.65 for 
all samples of Smo.55Sro.45Mn03 to simplify the problem 
somewhat. 

The normalized magnetization data for 
Smo.55Sro.45Mn03 |1S], has been scaled such that 
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FIG. 1: A graph of the Helmholtz free energy, for 5 = 0.65 
and X = 0.45, in the ferromagnetic and paramagnetic phases. 
The plot indicates tc « 0.2. 



is found to be A = 72.8nieV as 5 = 0.65. The value of 
Jpd found is in close agreement with the value of 60nieV 
as estimated for perovskite manganites by [50] . 



IV. SPECIFIC HEAT AND ENTROPY OF 

Smo.55Sro.45Mn03 

We can calculate the entropy of the system as: 
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In eq ( 10 1 Q is the total grand thermodynamic po- 
tential found from the CCDC theory. By substituting 
(7 = 1 or into equation! 10 1 and retaining the value of 



S as 0.65, we may calculate the entropy of the ferromag- 
netic and paramagnetic phases respectively. 

We suggest that the total entropy incorporating phase 
separation can be described as; 



the value of half volume fraction occurs at t = 0.2, see 
Fig(l2|. 
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FIG. 2: A graph of the volume fraction of Smo.55Sro.45Mn03, 
normalized such that the value of half volume fraction occurs 
at approximately, t = 0.2, as calculated for S = 0.65 and x = 
0.45. Solid circles represent scaled data [48], lines represent 
expressions used to describe the volume fraction. 

The transition temperature Tc in Smo.55Sro.45Mn03 
samples only varies by a few degrees Kelvin j48j , we can 
therefore look at the average conversion factor between 
the temperature in Kelvin and the normalized temper- 
ature t, which is 1.53 x 10^'^. This corresponds to an 
exchange energy between the carriers and localized spins 
of Jpd = 56meV, and hence the bipolaron binding energy 



Stot = VS/ + (l-V)Sp (11) 

where S / and Sp denote the entropy in the ferromagnetic 
and paramagnetic phases respectively, and V is the vol- 
ume fraction. 

Differentiating once more we may calculate the specific 
heat. 



We fit our calculated specific heat under constant vol- 
ume to the experimentally measured specific heat un- 
der constant pressure of Smo.55Sro.45Mn03, by assum- 
ing no thermal expansion. To fit to the experimen- 
tal data an estimate of the lattice specific heat is re- 
quired, typically the phonon blank method is used. The 
phonon blank method requires another sample to be pre- 
pared that does not exhibit the transition in question, 
this sample should have the same crystal structure and 
as close as possible a similar atomic mass. A possible 
phonon blank material for Smo.55Sro.45Mn03 could be 
produced by replacing Mn for Ti atoms, the resultant ma- 
terial would not exhibit the ferromagnetic-paramagnetic 
transition, if this material possessed the same crystal 
structure as Smo. 55810. 45Mn03 the specific heat of that 
sample would be very close to the lattice specific heat 
of Smo.55Sro.45Mn03. However no phonon blanks for 
Smo.55Sro.45Mn03 were produced. In fact we are un- 
aware of any published specific heat data for perovskite 
manganites including a phonon blank. We therefore sim- 
ply take an estimate of the smooth lattice background 
specific heat, see Fig(|3|. 

Figs([3]) and Q show that the theoretical specific heat 
is in reasonable agreement with the experimental results 
assuming our estimate of lattice specific heat is correct. 
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FIG. 3: Plots of specific heat of Smo.55Sro.45Mn03 annealed 
at (a)1200°C, (b)1300°C, (c)1400°C, (d)1500°C. Solid cir- 
cles represent experimental specific heat data |48l 149] . dashed 
lines correspond to our estimates of lattice specific heat, solid 
lines represent the addition of the lattice specific heat and the 
theoretical CCDC component. 



FIG. 4: Plots of the anomalous component of specific 
heat of Smo.55Sro.45Mn03 annealed at (a)1200°c, (b)1300°c, 
(c)1400°c, (d)1500°c. Solid circles represent the difference 
between the experimental specific heat data [481 149j and the 
estimate of the lattice specific heat, solid lines represent the 
specific heat calculated from the CCDC theory. 



We can calculate the entropy change associated with 
the anomalous part of the experimental data, FigQ, by 
means of; 




The entropy associated with the specific heat anomaly in 
the experimental data is plotted alongside the theoreti- 
cally determined entropy in Fig([5]). 

We investigate the entropy change from t = 0.0465 to 
t = 0.3, as these temperatures are well below and above 
the transition respectively, hence all samples should ex- 
hibit the same entropy change over this temperature in- 
terval. The only role that the phase separation plays is to 



broaden the observed specific heat peak. Table(|l]) shows 
the entropy change associated with the anomalous part 
of the data for each sample. 



Ta(°c) 


1200 


1300 


1400 


1500 


ASdata (JK-^mol-i) (l.dp) 


12.9 


12.8 


12.4 


13.0 



TABLE I: Entropy change associated with the anoma- 
lous component of the experimental specific heat data of 
Smo,55Sro.45Mn03, where Ta is the annealing temperature. 

Averaging the data from table(|l| we find the entropy 
change associated with the anomalous part of the exper- 
imental data from, t = 0.0465 to t = 0.3 as, ASdata = 
12.8 ± 0.1J/K~^mol~^. This is in good agreement with 
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V. RESISTIVITY OF Smo ssSm 45Mn03 

Using the same expressions governing the volume 
fraction, we fit the experimental resistivity data of 
Smo.55Sro.45Mn03 MM- 
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FIG. 5; Plots of entropy of Smo.55Sro.45Mn03 annealed at 
(a)1200°C, (b)1300°C, (c)1400°C, (d)1500°C. 



FIG. 6: Plots of the electrical resistivity of Smo.55Sro.45Mn03 
annealed at (a)1200°C, (b)1300°C, (c)1400°C, (d)1500°C. 
Solid circles show the experimental data |48l 149) ]. lines show 
fitted theoretical resistivity from equation! f4l. 



the theoretical entropy change over this temperature 
range, ASccdc = 12.802J/K"imor^ 

Several authors report that the entropy change as- 
sociated with the ferromagnetic-paramagnetic transition 
is much smaller than expected j49l I53j and hence un- 
explainable. Tanaka and Mitsuhashi'515 found a very 
small entropy change in Lao.8Cao.2Mn03, Ramirez 
states that this is approximately 10% of the expected 
value of « llJK~^mol~^ (for localized spins), approxi- 
mately twice as much entropy was found to reside under 
the ordering peak in Lao.eyCao.asMnOs by Ramirez et 
al i52j. The reasons for the entropy change being left 
unexplained up until now are that most authors [49l [53] 
have only considered local-moment degrees of freedom. 
The contribution to entropy of polarons acts to reduce 
the entropy change at Tc. Below Tc polarons are un- 
paired, above Tc the polarons are bound into immobile 
bipolarons, and hence the entropy of the charge carriers 
actually decreases at Tc. The total entropy change at Tc 
will always be positive, it is simply that the contribution 
of carriers acts to reduce this increase, resulting in en- 
tropy changes that are smaller than has been expected 



A simple analytical expression for the resistivity of a 
binary mixture, which is consistent with comprehensive 
numerical simulations was proposed by Kabanov et 
al |54j and considered for use, when explaining the 
resistivity of mixed phase perovskite manganites [43] : 



V „i-v 

r ferror para ' 



(14) 



which is valid in a wide range of the ratios 
Ppara/pferrolM], V is the volumc fraction. The resistivity 
is fitted in some- what of a generic manner, in as much 
that we are mixing the phases of a metal-like phase with 
that of an insulating phase. The expressions used are 
however, consistent with the physics of the CCDC model. 

Here we have used the expression governing the resis- 
tivity in the ferromagnetic phase as. 



Pfe 



aT*"" -I- c/sinh^ [fejo/keT] 



(15) 



In fitting to the experimental data the coefficient c is 
found as l(f2 cm) for all samples and hujo has been fixed 
as; 50meV, which is a typical value for the phonon cut 
off in LCMO \42i- 
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The expression for the ferromagnetic resistivity used 
here was suggested by Zhao et al[55]. Kubo and Ohata 
[55] have studied the magnon scattering for half met- 
als where the spin-up (minority) and spin-down (major- 
ity) bands are well separated. In this case, one-magnon 
scattering is forbidden, and thus the two-magnon pro- 
cess plays a role in the low-temperature resistivity, which 
gives a contribution that is proportional to T"*'^. The 
term c/sinh^ [fiw/ksT] comes from the relaxation time 
of polarons calculated by Lang and Firsov [57]. Jamie 
et al [55] have ruled out electron-electron, and electron- 
phonon scattering mechanisms in the low temperature 
conductivity of manganites, that would give T^ and 
contributions respectively. 

The resistivity of the paramagnetic phase can be well 
described by [43) : 

Ppara = bTe^ (16) 

We also include a residual resistivity po- The parame- 
ters used to fit the resistivity are shown in table ([TT|). 



Ta (°C) 


1200 


1300 


1400 


1500 


pa (ri cm) 


3 


1 


0.21 


0.2 


a (n cm K^"'^) X 10""' 


38 


17 


4 


0.07 


b {Q. cm K-i) X 10-*^ 


2.2 


1.3 


1.25 


1.2 


E (meV) 


138 


131 


123 


116 


Ea (meV) 


101.6 


94.6 


86.6 


79.6 



TABLE II: Parameters governing the resistivity of 
Smo.55Sro.45Mn03 l48ll49j 

It is clear that the peaks observed in the resistivity of 
Smo.55Sro.45Mn03, FigQ, are broadened with increasing 
disorder and phase separation. One can also see that the 
resistivity is increased with the increasing disorder, this 
is simply due to increased scattering rates, as can be 
seen by the increase values of the constants a,b and po in 
table(|ll]) with decreasing annealing time. 

The high temperature mobility of polaronic carriers is 
likely to be dominated by hopping events due to polaron 
band narrowing, [571 159) , therefore E has contributions 
from the bipolaron binding energy, A, and the activation 
energy of polarons, Ea as; 

E = Ea + A/2 (17) 



By fitting the specific heat we have found A/2 = 
36.4meV, and from fitting the resistivity we have found 
E. This has allowed us to separate the polaron hopping 
activation energy and bipolaron binding energy from the 
exponential behaviour of the resistivity in the paramag- 
netic regime. See table([ll]) for values of polaron hopping 
activation energy. We can see that the hopping activation 
increases with an increasing density of grain boundaries 
(decreasing annealing temperature) [H] . This is what one 
would expect to see as the increasing levels of disorder 
act to localize the carriers resulting in a larger activation 
energy required for hopping. 

VI. CONCLUDING REMARKS 

Introducing Phase separation into the thermodynamics 
of the microscopic CCDC model by mixing the high and 
low temperature phases has allowed for the specific heat 
of Smo.55Sro.45Mn03 samples produced by Egilmez et al 
[48| |49] to be well explained. Up until now the entropy 
change associated with the ferromagnetic-paramagnetic 
transition in colossal magnetoresistive perovskite man- 
ganites has been left largely unexplained, however we find 
the entropy change from the CCDC theory and the exper- 
imental data to be in close agreement. The contribution 
to entropy of polarons acts to reduce the entropy change 
at Tc. Below Tc polarons are unpaired, above Tc the po- 
larons are bound into immobile bipolarons, and hence the 
entropy of the charge carriers actually decreases at Tc. 
The total entropy change at Tc will always be positive, it 
is simply that the contribution of carriers acts to reduce 
this increase, which results in an entropy change that is 
smaller than has been expected [121 [53] . Fitting specific 
heat data, in conjunction with a generic phase separated 
model of resistivity has allowed bipolaron binding energy 
and polaron hopping activation energy contributions to 
the exponential behaviour of resistivity to be separated. 
We hope that experimentalists working in the field of 
perovskite manganites try to produce phonon-blank ma- 
terials such that the lattice contribution to specific heat 
of samples can be accurately determined. 

We accept that the CCDC theory is not the only the- 
ory that hopes to explain the physics of colossal magne- 
toresistive materials, we have shown however for the first 
time that this theory may be used to explain the specific 
heat of perovskite manganites. More work is required, es- 
pecially by experimentalists to either conclusively prove 
or disprove this model. 
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